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Abstract 

The nematic-to-isotropic orientational phase transition, or equivalently the RP 2 model, is consid- 
ered in two dimensions and the question of the nature of the phase transition is addressed. Using 
powerful conformal techniques adapted to the investigation of critical properties of two-dimensional 
scale-invariant systems, we report strong evidences for a transition governed by topological defects 
analogous to the Berezinskii-Kosterlitz-Thouless transition in two-dimensional XY model. 



Liquid crystals may be seen as constituted of 
molecules essentially represented by long rigid rods. 
From maximization of entropy at high tempera- 
tures, all the molecule orientations are equally prob- 
able, independently of the neighbouring molecule 
directions and the system exists in an isotropic 
phase. At low temperatures a preferential orien- 
tation is more favourable in order to minimize in- 
teraction terms, and an ordered structure emerges. 
When order occurs along one space dimension only, 
the system is said to be nematic. Still at lower 
temperatures, other ordered phases can appear, e.g. 
smectic phases. 

In a lattice model, each molecule may be repre- 
sented by a unit vector a w at site w of an hyper- 
cubic lattice A of linear extent L. The <r's live in 
a three-dimensional space attached to each lattice 
site. In the nematic phase, the preferential direc- 
tion defines a unit vector, n, called the director, 
and one can measure the deviation of molecule <x w 
with respect to the director by the scalar product 
a w ■ n = cos 0m . Due to the local Z 2 symmetry 
(the rods are not oriented), one cannot distinguish 
between opposite directions 9 W and 9 W + n, and 
cos 9 W vanishes on average while cos 2 9 W does not. 



In the disordered phase on the other hand, the an- 
gles are measured with respect to any arbitrary di- 
rection, and the thermal average of course leads to 



(cos 9 W ) = and (cos 2 9 V 



3, so that (cos 2 9 W ) — | 
represents a convenient order parameter. In the lit- 
erature on liquid crystals, one usually defines the 
local order parameter by the second Legendre poly- 
nomial, 



m[w) 



(P 2 (* w -n)) = (P 2 {cos9 w )}. (1) 



This definition suggests to consider the following 
Hamiltonian to describe the nematic transition, 



H 



J 



k R T k R T 



(2) 



where [i stands for the unit basis vectors of the lat- 
tice, cr w ■ (T w+il = cos(9 w — 9 W+Ii ) is the scalar prod- 
uct between neighbouring vectors distant from one 
lattice spacing, and the interaction term —JPi{(j w ■ 
cr w+fl ) is reminiscent from a dipole-dipole interac- 
tion. This Hamiltonian was introduced by Leb- 
wohl and Lasher [1] lattice version of the mean 
field theory of Maier and Saupe [2], and its success 
came from its ability to reproduce the weak first or- 
der phase transition observed experimentally in the 
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three dimensional nematic transition [3]. In a more 
abstract context, this Hamiltonian is known as the 
RP 2 model, since at each lattice site is attached the 
manifold of directions in 3-dimensional space, also 
called the real projective space in 3 dimensions [4]. 

Like the non linear cr-model, this model pos- 
sesses generically the symmetry group 0(n) which is 
non abelian for n > 3, and specifically the Lebwohl- 
Lasher or RP 2 model has a 0(3) symmetry. The 
question, as it was outstandingly formulated by 
Kunz and Zumbach [4], of the nature of the tran- 
sition in the RP 2 model in two dimensions is still 
incompletely solved. The existence of a phase tran- 
sition (at finite temperature) in two dimensional 
systems seems connected to the abelian nature of 
the underlying symmetry group, as both Ising and 
XY models are famous examples, unlike the Heisen- 
berg model. On the other hand, according to the 
Hohenberg-Mermin- Wagner theorem [5, 6], models 
possessing a continuous symmetry group cannot ex- 
hibit any finite macroscopic magnetization with no 
magnetic field applied in dimensions 1 or 2 (we in- 
tentionally use the familiar terminology of magnetic 
systems). The two-dimensional XY model is the 
most famous example and it exhibits a non con- 
ventional transition [7, 8, 9]. In spite of the ab- 
sence of long-range order (LRO) at low tempera- 
tures, the spin-spin correlation function decays al- 
gebraically with an exponent which increases mono- 
tonically with temperature [10, 11, 12] up to a 
temperature called after Berezinskii, Kosterlitz and 
Thouless [13, 14, 15]. In this critical phase, macro- 
scopic ordering is prevented by collective excita- 
tions, namely spin waves which nevertheless do not 
exclude a coherent orientation of spins at a smaller 
length scale. Together with the spin waves, local- 
ized excitations appear with increasing tempera- 
ture. These are topological defects associated in 
pairs, like pairs of opposite charges in the low tem- 
perature phase of a two-dimensional Coulomb gas, 
and they perturb the spin field only locally. This 
is the usual meaning of the term quasi-long-range 
order (QLRO). The Berezinskii-Kosterlitz-Thouless 
transition is governed by unbinding of these defects 
which completely suppresses any type of long-range 
order at high temperature, hence the correlation 
functions decay exponentially like in an ordinary 



paramagnetic phase. This is the standard scenario 
of a topological transition. The puzzle becomes con- 
fused when one notices that XY or Heisenberg mod- 
els in three dimensions display conventional con- 
tinuous transitions and the RP 2 model exhibits a 
first order transition while renormalization group 
treatment of non-linear a model predicts the ab- 
sence of any transition at non-zero temperature in 
2d (asymptotic freedom in the context of lattice 
gauge theories [16]) and a continuous one in 3d for 
all the three models. We also have to mention that 
the RP 1 model (the same as given in equation (2), 
but with two-component vectors cr w ) exactly coin- 
cides with the XY model. The question of a pos- 
sible topological transition in the two-dimensional 
(non abelian) Lebwohl-Lasher model is thus par- 
ticularly attracting and was already addressed in 
previous studies [17, 18, 19, 20, 4]. In their re- 
markable work, Kunz and Zumbach [4] concluded 
in 1992 in favor of such a topological transition sce- 
nario, essentially on the basis of qualitative argu- 
ments (pairing of topological defects at low tem- 
perature where they carry most of the energy in 
the system, sharp increase of the density of defects 
and apparent discontinuity of the rotational rigid- 
ity modulus at the transition, finite cusp in the spe- 
cific heat, proliferation of unbinded defects at high 
temperature). Even though they performed a care- 
ful and sizeable study, they were unfortunately not 
able to decide conclusively between essential singu- 
larities or standard power laws - though their pref- 
erence was for the first case - for the correlation 
length and the susceptibility when approaching the 
transition from the high temperature phase. This 
is essentially due to the limited possibilities of com- 
puters ten years ago, since the authors already took 
care about potential critical slowing down problems 
as they adapted the Wolff cluster algorithm to the 
R pn-i mo del. 

Ten years later, we want to address the same 
question of the nature of the phase transition of the 
two-dimensional Lebwohl-Lasher model. Since we 
believe that the conclusions of Kunz and Zumbach 
will hardly be improved, even with more powerful 
facilities, it is necessary to reconsider the problem 
from a different point of view. We will thus assume 
the existence of a critical phase at low tempera- 
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tures and then follow the line of the behaviour of 
the XY model to predict consequences of the above 
mentioned assumption, consequences which may be 
compared easily to numerical results. Of course, the 
test must discriminate between different scenarios, 
starting from the confirmation (and thus the char- 
acterization) of a topological transition, a conven- 
tional continuous transition, a first-order one, or no 
transition at all. The existence of a scale-invariant 
low temperature critical phase is characteristic from 
the first situation. Such a system must thus be con- 
formally invariant at any temperature below the 
transition Tkt (we will abusively keep the termi- 
nology adapted to the XY model), so it becomes 
advantageous to deduce the functional expression 
of the correlation functions or density profiles in a 
restricted geometry adapted to numerical simula- 
tions from a conformal mapping w(z): 

(<t w1 -<t W2 ) = \v/(z 1 )\~ x '\v/(z 2 )\~ x '{(r gl -(T ga ). (3) 

Here, w labels the lattice sites in the transformed 
geometry (the one where the computations are re- 
ally performed), z are the corresponding points in 
the original one (usually the infinite plane where 
} ~ \zi — z 2 \ VtT takes the standard power- 
law expression), and x a = \r] u is the scaling dimen- 
sion associated to the scaling field a. The interest 
of such an approach lies in the full inclusion in the 
functional expression of the changes due to shape 
effects. The most famous example is the exponen- 
tial decay of the correlation functions at criticality 
along a strip of finite width, unlike the algebraic de- 
cay in the infinite plane. For simplicity reasons, it is 
even more convenient to work with density profiles 
in a finite system with symmetry breaking fields 
along some surfaces in order to induce a non van- 
ishing local order parameter in the bulk. In the 
case of a square lattice A of size L x L, with fixed 
boundary conditions along the four edges <9A, one 
expects [21, 22, 23] 

m FBC (w) = (P 2 (Cit, • h 9A(u ,))) F Bc ~ [k(u>)]"%) 



x |(l- S n 2 2K2£) (1 - fc 2 sn 2 2K2£)| 2(5) 

where Fbc specifies that the boundary conditions 
are fixed. This expression easily follows from the 



expression of the order parameter profile decaying 
in the upper half-plane from a distant surface of 
spins constantly fixed in a given direction, m{z) = 
{P2{tr z • haA.( 2 )))uHP ~ y~ Xa ', and from the confor- 
mal transformation of the upper half-plane (uhp) 
z = x + iy(0<y< 00) inside a square w = u + iv 
of size L x L (— L/2 < u < L/2, < v < L) with 
open boundary conditions along the four edges, re- 
alized by a Schwarz-Christoffel transformation [24] 

W ^ = 2K z = sn \ir)' ^ 

Here, F(z, k) is the elliptic integral of the first kind, 
sn(2Kw/L) the Jacobian elliptic sine, K = K(k) = 
F(l,k) the complete elliptic integral of the first 
kind, and the modulus k depends on the aspect ra- 
tio of A. 

Our strategy is now to fit numerical data of the 
order parameter profile against expression (4). Like 
in the previous study of Kunz and Zumbach [4, 20] 
the resort to a cluster update algorithm is neces- 
sary in order to prevent the critical slowing down, 
all the spins of clusters (build through intermedi- 
ate bond variables) being updated simultaneously. 
The algorithm becomes particularly efficient if the 
percolation threshold of the bond process occurs 
at the transition temperature of the spin model, 
which ensures the updating of clusters of all sizes 
in a single MC sweep. For 0{n) models, Ising vari- 
ables are defined in the Wolff algorithm by the sign 
of the projection of the spin variables along some 
random direction. The bonds are then introduced 
through the Kasteleyn-Fortuin random graph repre- 
sentation [25]. When one uses fixed boundary con- 
ditions, a difficulty occurs and the Wolff algorithm 
should become less efficient, since close to criticality 
the unique cluster will often reach the boundary and 
no update is made in this case. This is circumvented 
by the following trick: even when the cluster reaches 
the fixed boundaries, it is updated - and so are the 
boundary spins - and the order parameter profile is 
then measured with respect to the new direction of 
the boundary spins, m FBC (w) = (P 2 (<t iu • <t 9A )) Fbc . 

Using this procedure, we studied systems of size 
48 x 48 up to 200 x 200. For the measurement of 
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the order parameter profile, we discarded 10 6 Wolff 
sweeps for thermalization, and the measurements 
were performed on 10 6 production sweeps. For rea- 
sons which are made obvious below, the energy den- 
sity required a better statistics and the measure- 
ments were obtained over 16. 10 6 sweeps. 

In order to underscore the existence of a line of 
marginal fixed points in the low temperature phase, 
we first check qualitatively the expression of the 
energy-energy correlations from the behaviour of 
the energy density profile. The energy density at 
site w is for example defined as the average value of 
the energies of the four links reaching w: 



(7) 

The existence of a regular contribution in the energy 
density makes the calculation a bit more subtle than 
what presented in equation (4). This regular con- 
tribution {eq(T)) which depends on T cancels after 
a suitable difference between profiles obtained with 
different conditions at the boundaries (free (fBC) 
and fixed (Fbc) boundary conditions). Although 
it makes the numerical computation longer in or- 
der to reach some satisfactory accuracy, this makes 
possible to extract the singularity associated to the 
energy density: 



(£z)fbc 

(£z)fBC 



(e (T)) + B FBC (T)y-^ T )/ 2 , (8) 
(e (T)) + B iBC (T)y-^ T )/2. (9) 



This is clearly illustrated in figure 1 where con- 
vergence towards the same temperature-dependent 
constant {£q(T)) is shown, with amplitudes of the 
singular terms having opposite signs, therefore a 
simple difference of the quantities measured in the 
square geometry, 

Ae(w) = (e w ) FBC - (e w ) {BC ~ AB x [k(w)] ~M t ) 

(10) 

leads to the value of the thermal scaling dimension 
r/ £ (T). The right part in figure 1 presents a log-log 
plot of the difference Ae(w) vs k(w) at two temper- 
atures below Tkt and one above which shows that 
the functional expression used is no longer valid, as 
expected, in the paramagnetic phase. Due to the 
strong fluctuations, in the QLRO phase the data 



scatter around straight lines which represent the 
slopes {k(w)]~ 2 . This figure, though not definitely 
conclusive, confirms that the exponent of the de- 
cay of energy-energy correlations keeps a constant 
value r] £ {T) = 4 in the low-temperature phase of 
the RP 2 model, confirming that like in the case 
of the XY model, the temperature is a marginal 
field, responsible for the existence of a critical line 
in the whole low-temperature phase. It thus implies 
a thermal scaling exponent x £ = d—yt = 2 which en- 
sures a vanishing RG eigenvalue yt = (up to Tkt 
where it is consistent with an essential singularity 
of £ above the KT point, as suspected by Kunz and 
Zumbach [4]). The energy-energy correlation func- 
tion in the plane should thus decay algebraically as 



{e Zl e Z2 ) ~ \zi - z 2 \ 
with r? e (T) = 2x £ = 4 VT < Tkt- 



11 



i 

0.95 
0.9 



A 0.85 
V 0.8 



7=0.1- 
T=0.2- 

T=0.3- 



T=0.6 





0.1 0.2 0.3 0.4 0.5 




v,, T=0.6 



N 



T=0A tta^J 



.-4< 



k(w) 



k(w) 



Fig. 1: MC simulations of the 2c? RP 2 model inside 
a square of 100 x 100 spins (16. 10 6 MC sweeps after 
cancellation of 10 6 for thermalization). Several temper- 
atures below the transition temperature are shown, and 
one above (T = 0.6). Left: local energy density vs the 
rescaled variable k(w). Right: log-log plot of the differ- 
ence Ae(w). 

From the existence of a marginal line, one may 
suspect that the other scaling dimensions should 
continuously vary with the marginal field. The or- 
der parameter profile thus has to obey equation (4) 
in the whole low temperature phase, but with an 
exponent T] a (T) which depends on T. Equivalently, 
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a log-log plot of mpBc(w) vs k(w) must display 
straight lines with different slopes below Tkt- This 
is exactly what is observed in figure 2. Again, the 
curves start to deviate from the straight line when 
the system enters the high temperature phase. 
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Fig. 2: MC simulations of the 2d RP 2 model inside a 
square of 100 x 100 spins (4.10 6 MC sweeps after cancel- 
lation of 10 6 for thermalization) . The order parameter 
density is plotted against the rescaled distance at differ- 
ent temperatures below the KT transition temperature 
and slightly above. The numbers on the right give the 
value of ksT/J, and 0.537 (the estimate of Kunz and 
Zumbach [4]) appears to be overestimated. 
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Fig. 3: Left: Correlation function decay exponent of 
the 2d RP 2 model as a function of the temperature. The 
vertical dashed line is the estimate of the transition tem- 
perature of Kunz and Zumbach [4]. Right: value of the 



X 2 per degree of freedom resulting from the power law 
fits of the data in figure 2. The sharp increase is an 
indication of the location of the transition temperature 
(already below the estimate of Kunz and Zumbach [4]). 

The slopes measured in figure 2 are reported as 
a function of the temperature in figure 3. 
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Fig. 4: Dependence of the rotational rigidity modulus 
with the rj a exponent. 

Another signature of this mechanism is the be- 
haviour of the rotational rigidity modulus. This 
latter quantity generalizes the helicity modulus T 
in the XY model, which measures the quadratic re- 
sponse in the free energy of the system to a twist 
accross the sample. This is generalized to the RP 2 
model by measuring the change in free energy when 
a rotation of angle eft around some axis in spin 
space is applied to the system: F{4>) — F(Q) = 
T,cp 2 -\-0((j) A ). This expression defines the rotational 
rigidity modulus S. In the XY model, there ex- 
ists a universal relation between the helicity mod- 
ulus and the correlation function exponent, r/ a = 
ksT '/2ttJT [26], which appears as a consequence of 
the Kosterlitz recursion relations. In the Lebwohl- 
Lasher model, the same type of behaviour is checked 
in figure 4 where one observes a linear dependence 
(the larger the system size, the better the linear 
behaviour) of r/ CT with T/E, correctly fitted at low 
temperatures by i] a = 0.117 fc^T/JS. 

To conclude, we mention that the most remark- 
able feature is a complete analogy with the two- 
dimensional XY model [22] where the transition is 
mediated by the defects. The exponent r\ e keeps a 
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constant value while rj a , associated to the order pa- 
rameter, starts from zero at T = and increases 
linearly with T at low temperatures where a spin 
wave approximation should capture the essentials 
of the behaviour of the system. The influence of 
pairs of topological defects, which would appear in 
increasing number when the temperature increases, 
is probably responsible for the deviation from the 
spin wave approximation and of the sharper in- 
crease of r] a , and the transition is presumably com- 
pletely governed by unbinding of these defects, like 
in the Kosterlitz-Thouless scenario. The order pa- 
rameter exponent at the transition /cbTkt/-^ — 0.52 
takes a value r) a (T) ~ 0.40(2). The relation between 
rotational rigidity modulus and the exponent rj a (T) 
also seems completely coherent with what happens 
in the XY model, giving one more evidence of the 
topological nature of the transition. 
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